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B+(1S)  +  H2  ==>  BH+(2Z)  +  H 

A  Woodward-Hoffmann  Forbidden  Ion-Molecule  Reaction 


Jeff  Nichols#,  Samuel  J.  Cole+,  Maclej  Gutowski,  and  Jack  Simons 
Chemistry  Department 
University  of  Utah 
Salt  Lake  City,  Utah  84112 

Abstract 

The  reaction  of  B+^S)  with  H2  on  the  ground  potential  energy  surface  is 
examined  using  ab  initio  electronic  structure  methods.  A  weakly  bound  T-shaped 
B+— H2  complex  of  C2v  symmetry  is  found  to  lie  855  cm'1  below  the  B+  +  H2  reactant 
energy.  Its  H-H  internuclear  distance  is  only  slightly  distorted  from  that  of  H2;  the  B-H 

distance  (ca.  2.6  A)  is  much  longer  than  the  covalent  bond  length  in  BH+  (1.2  A). 

Further  along  the  reaction  coordinate  is  found  a  narrow  'entrance  valley'  characterized 
by  strong  B+-to-H2  inter-reactant  forces  but  very  small  distortion  of  the  H-H  bond 
length  or  the  H-H  vibrational  frequency.  As  one  proceeds  further  up  this  valley,  a 
geometry  is  reached  at  which  the  asymmetric  distortion  mode  (of  b2  symmetry) 
develops  negative  curvature.  Distortion  along  the  asymmetric  mode  leads  to  a  transition 
state  of  Cs  symmetry  posessing  one  'long'  B-H  distance  (rBH=  1-765  A),  one  'short'  B- 
H  distance  (reH'  =  1 .251  A)  and  an  essentially  'broken'  H-H  bond  (rHH  -  1-516  A).  Its 
energy  lies  22,528  cm'1  above  B+  +  H2  and  2031  cm*1  or  ca.  0.25  eV  above  the 
thermodynamic  reaction  threshold  for  BH+  +  H  formation,  which  is  predicted  to  be 
endoergic  by  20,497cm*1  •  A  geometrically  stable  linear  HBH+(1Z)  species  is  found  to 
lie  19,259  cm*1  below  B+  +  H2.  Its  BH  bond  length  (rBH=1-173  A)  is  only  slightly 
longer  than  that  in  BH+  (1 .2  A).  Similarities  and  differences  among  reactions  of  the 
isoeiectronic  Be,  B+,  and  Li*  with  H2  are  also  discussed. 
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I.  Introduction 


Earlier  theoretical  calculations1  on  Be(1S)  +  H2  ==>  HBeH  (1l)  and  on2 
Mg(1S)  +  H2  =«>  HMgH(1z)  yielded  qualitatively  similar  ground-state,  C2v  - 
constrained  potential  energy  surfaces.  A  two-dimensional  contour  characterization  of 
such  C2v  -  constrained  potential  energy  surfaces  is  given  in  Fig.  1.  Both  of  these 

alkaline  earth  atoms,  as  well  as  the  analogous  B+  ion  considered  here,  have  ns2  1 S 
ground  electronic  states  and  relatively  low-lying  nsnp  1-3P  excited  states.  As  discussed 
in  detail  below,  the  chemical  reactivity  of  the  1 S  ground-states  of  such  species  is 
strongly  influenced  by  the  presence  of  the  excited  1  >3P  states.  Before  discussing  our 
motivations  for  undertaking  the  theoretical  study  of  the  B+(1S)  +  H2  ==>  BH+(2£)  + 

H  reaction  and  our  findings,  let  us  exploit  our  experience  on  the  Be  and  Mg  reactions  to 
anticipate  some  of  the  electronic  structure  characteristics  that  would  be  expected  to 
affect  strongly  the  dynamics3a>3b.3c  (e.g.,  energy  threshold,  cross-section  magnitude, 
and  energy  dependence)  of  the  B+  +  H2  reaction. 

A.  The  Electronic  Configurations  of  Interest  in  C2v  Symmetry 

1.  The  M  +  H2  Reactant  and  HMH  Product  Configurations 

Briefly,  the  ns2  configuration  of  the  reactant  M  atom  (or  ion)  when  combined 
with  the  a2g  configuration  of  the  H2  molecule  in  its 

X  1lg  ground  state,  does  not  correlate  with  the  ground-state  o2g  o2u  configuration  of  the 
linear  HMH  molecule  (or  ion)  In  C2v  symmetry,  the  M  +  H2  reactants  have  the  ns2  a2 g 
=■>  1a2-|2a2i  electronic  configuration,  and  the  HMH  products  have  the  o2g  o2u  = 
1a2-|1b22  configuration.  Moreover,  neither  of  these  two  configurations  can  describe  the 
MH  +  H  products  which  require  an  'open-shell'  configuration  of  the  form  o2mh  om  oh- 
For  the  M  +  H2  reactants,  the  labels  lai,  2ai,  and  1b2  are  used  to  denote  the  valence  M 
ns,  H2  og,  and  M  (in  plane)  np  orbitals.  The  M  inner-shell  orbitals  are  not  explicitly 
identified  in  this  abbreviated  notation.  For  the  HMH  product  molecule,  1  a-| ,  2ai ,  and 
1b2  denote  the  M-H  o,  o*.  and  o  orbitals,  respectively.  For  the  MH  +  H  products,  omh 
labels  the  M-H  bonding  0  orbital,  om  the  M-centered  non-bonding  orbital,  and  oh  the 
H-atom  Is  orbital. 

Both  of  the  above  dominant  configurations  of  the  M  +  H2  reactant  and  MH2 
product  are  of 1  Ai  symmetry,  so  they  mix  as  one  proceeds  along  the  reaction  path  to 
produce  the  'avoided  crossing'  which  gives  rise  to  the  well  known  symmetry  imposed 
activation  barrier  characteristic  of  such  Woodward- Hoffmann  forbidden  reactions.  For  a 
C2v  symmetry  preserving  reaction  path  appropriate  to  the  insertion  of  B+  into  H2  to 
produce  HBH+,  these  correlations  are  described  semi-quantitavely  in  Fig.  2. 


2.  Low-Energy  Excited  State  Configurations 


The  excited  1  »3P  states  of  the  M  species,  when  interacting  with  H2  in  C2v 
symmetry,  give  rise  to  singlet  and  triplet  states  of  A-|,  Bi,  and  B2  symmetry.  The 
1  >3B2  states  posesses  the  most  attractive  interactions  because  they  allow  M's  in-plane 
np  orbital  of  b2  symmetry  to  interact  constructively  with  the  H2  molecule's  'empty' 
antibonding  ou  orbital  as  shown  in  Fig.  3.  The  relative  locations  of  these  1  *3B2  states 
along  the  insertion  reaction's  C2v  reaction  path  are  also  depicted  in  Fig.  2  (with  the 
triplet  states  lying  below  the  singlets). 

B.  Implications  for  Reactivity 

The  fact  that  the  ‘•Ai  ground  state  surface  is  either  intersected  or  closely 
approached  by  the  1  -3B2  surfaces  has  important  consequences  for  the  insertion  reaction 
under  discussion.  In  particular,  collsions  entering  the  region  where  the  1  Ai  and  1 B2 
surfaces  are  close  in  energy  may  'hop'  from  the  ns2  1  Ai  entrance-channel  surface  to  the 
1 B2  surface  if  the  collision  occurs  with  some  degree  of  assymmetry.  Collisions  that 
occur  slightly  away  from  C2v  symmetry  will  experience  weak  coupling  of  the  and 
1B2  surfaces  (which  are  both  of  1A'  symmetry  in  the  lower  Cs  point  group).  The  point 
is  that  any  loss  of  C2v  symmetry  permits  the  ^Ai  and  1B2  states  to  mix,  thereby 
allowing  reactants  to  move  onto  the  1 B2  surface,  which,  in  Cs  symmetry,  correlates  to 
MH  (2I)  -1-  H  products.  Of  course,  if  any  appreciable  spin-orbit  coupling  is  operative, 
transitions  to  the  3B2  surface,  which  also  correlates  directly  to  the  MH  (2z)  +  H 
products,  can  also  occur.  It  is  through  such  surface  'hoppings'  that  the  MH  +  H  channel  is 
accessed  when  the  reaction  begins  with  ground-state  M  and  H2  species. 

C.  Reaction  Path  "Shape' 

A  somewhat  more  quantitative  view  of  the  C2v  ground-state  energy  surface 
appropriate  to  these  'insertion'  reactions  is  provided  by  Fig.  1  for  the  ab  initio 
calculated1  1  Ai  surface  of  Be(1S)  +  H2  ==>  HBeH  (X  1lg).  This  surface  and  others  of 
this  M  (ns2;  1S)  +  H2  ==>  HMH  ^Xg)  family  are  characterized  by  potential  energy 
landscapes  along  which: 

a.  The  M  +  H2  entrance  channel  is  very  'straight'  (i.e.,  the  reaction  coordinate  is 
dominated  by  M-to-H2  relative  motion  with  very  little  H-H  displacement)  and  has  very 
large  positive  curvature  transverse  to  the  reaction  coordinate  (i.e.,  the  H-H  bond 
remains  intact  and  is  very  ’stiff). 

b.  Once  a  critical  M-to-H2  distance  is  reached,  the  reaction  coordinate  undergoes  a 
sudden  change  to  become  dominated  by  H-H  stretching  with  much  less  M-to-H2 
movement. 

c.  Further  along  the  reaction  path,  a  C2v  transition  state  is  reached  that  lies  more  than  2 
eV  above1  the  energy  of  the  M  +  H2  reactants  and  even  above  the  energy  of  the  MH  +  H 
fragments  that  could  be  formed  if  C2v  symmetry  were  not  (artificially)  enforced. 

d.  Between  the  point  where  the  reaction  path  acquires  negative  curvature  along  its 
'uphill'  coordinate  and  the  C2v  -  constrained  transition  state,  a  region  appears  within 
which  negative  curvature  also  exists  along  the  asymmetric  (b2)  distortion.  Because  of 
the  constraint  to  C2v  symmetry,  the  force  (i.e.,  teh  energy  gradient)  along  this  b2 
direction  vanishes  identically.  The  appearance  of  negative  curvature  along  a  direction 
transverse  to  the  reaction  coordinate  indicates  that  a  lower-energy  path  can  be  found  if 


one  moves  away  from  the  C2v  geometry  along  this  asymmetric  direction.  In  the  above 
cases,  doing  so  eventually  leads  to  formation  of  ground-state  MH  +  H  products.  The 
negative  curvature  on  the  ‘•Ai  surface  along  the  b2  direction  is  a  signal  that  a  ^2 
surface  is  in  close  proximity;  it  does  not  mean  that  an  intersection  with  a  1 B2  surface  is 
taking  place,  but  that  such  a  surface  is  energetically  nearby.  It  is  in  this  region  of  the 
surface  that  one  must  break  C2v  symmetry  in  search  of  the  true,  lowest-energy, 
transition  state. 

e.  If  C2v  symmetry  constraints  are  kept  operative,  movement  along  the  b2  asymmetric 
mode  will  not  occur,  and  the  reaction  coordinate  will  evolve  smoothly  to  and  beyond  the 
C2v  -  constrained  transition  state  mentioned  above. 

f.  Beyond  this  C2v  transition  state,  the  reaction  coordinate  eventually  develops  positive 
curvature  as  the  linear  HMH(1lg)  geometry  is  approached. 

D.  Why  Study  B+C  S)  +  H2? 

The  present  work  is  aimed  at  extending  the  investigations  described  above  in  at 
least  two  aspects:  (i)  to  include  a  positive  metal  ion  (but  still  with  a  closed-shell  1 S 
ground  electronic  state)  as  the  reactive  species  so  any  effects  caused  by  ion-molecule 
interactions  can  be  examined,  and  (ii)  to  investigate  reaction  paths  that  are  not  C2v  - 
preserving  so  as  to  permit  the  MH  +  H  product  channel  to  be  explored.  We  chose  B+^S) 
as  the  positive  ion  because:  (i)  it  is  isoelectronic  with  Be(1  S)  which  we  examined 
earlier  (and  thus  similarities  and  differences  between  the  Be  and  B+  cases  are  of 
interest),  (ii)  it  allows  highly  accurate  calculations  to  be  performed  with  reasonable 
effort,  and  (iii)  experimental  guided  ion  beam  and  other  data3a-3d  giving  the  cross- 
section  for  BH+  production  as  a  function  of  B+  kinetic  energy  are  available  and  in  need 
of  theoretical  interpretation. 

II.  Computational  Methods 

A.  Basis  Sets 

The  basis  set  for  the  H  atoms  consists  of  the  Dunning  augmented  correlation 
consistent  (cc)  polarized  valence  triple-zeta  (p-VTZ)  [5s2p1d|  3s2p1d]  set4  of 
functions.  For  the  B+  ion,  the  Dunning  [10s5p2d|  4s3p2d]  augmented  cc  p-VTZ  basis 
set  was  used.  A  total  of  55  contracted  Gaussian-type  basis  functions  resulted. 

B.  Electronic  Configurations  and  Wavefunctions 

Both  multiconfigurational  self-consistent  field  (MCSCF)  and  coupled-cluster 
methods  were  used  to  treat  correlations  among  the  four  valence  electrons  of  the  BHH+ 
system.  In  particular,  the  CCSD(T)  variant5  of  the  coupled-cluster  approach6,  which 
includes  all  single  and  double  excitations  in  a  fully  correct  manner  and  treats  triple 
excitations  by  approximate  non-iterative  means,  was  employed. 

The  discussion  of  Sec.  I  makes  it  clear  that  no  single  electronic  configuration  can 
describe  even  the  ground  state  of  this  system  throughout  the  C2v  or  C$  reaction  paths. 

For  this  reason,  multiconfigurational  methods  were  required.  In  the  MCSCF  calculations, 
the  four  valence  electrons  were  distributed,  in  all  ways  consistent  with  overall  spatial 
and  spin  symmetry,  among  6  valence  orbitals.  In  the  C2v  calculations,  4  of  these 
orbitals  were  of  at  symmetry,  1  was  of  b2  symmetry,  and  1  was  of  bi  symmetry.  In  C$ 
symmetry,  there  were  5  of  a’  symmetry  and  1  of  a"  symmetry.  The  two  B+  Is  electrons 


were  constrained  to  occupy  a  single  a-|  or  a'  core  orbital  in  all  of  the  electronic 
configurations  generated.  This  process  generated  41  electronic  configurations  of 1  A-| 
symmetry  in  the  C2 v  point  group  and  65  of 1  A'  symmetry  in  the  Cs  point  group. 

The  above  MCSCF  calculations  were  employed,  along  with  our  Utah  MESSKit7 
analytical  energy  derivative  and  potential  energy  surface  'walking'  algorithms8  to  find 
and  characterize  (via  geometry  and  local  harmonic  vibrational  frequencies)  the  local 
minima,  transition  states,  and  reaction  paths  discussed  below.  Near  each  such  point,  the 
CCSD(T)  method  was  used  to  evaluate  the  total  correlation  energies  at  a  more  accurate 
level;  finite  difference  methods  were  also  employed  within  the  CCSD(T)  approach  to 
refine  the  predicted  geometry  of  each  such  point  on  the  surface. 

III.  Findings  and  Comparison  to  Guided  Ion  Beam  Results 

As  detailed  in  Table  I  and  described  qualitatively  in  Figs.  2  and  4,  we  find  the 
potential  energy  surface  for  B+^S)  +  H2  ==>  BH+(2I)  +  H  to  possess  most  of  the 
features  that  are  expected  from  the  discussion  in  Sec.  I.  Although  there  are  differences 
that  might  have  been  expected  because  of  the  additional  long-range  ion-molecule 
interactions  that  are  operative  in  this  case,  we  find  such  affects  to  be  quite  small  in  this 
case.  Our  primary  findings  are  summarized  as  follows: 

a.  A  weakly  bound  B+— H2  complex  lies  1143  cm*1  below  the  B+  +  H2  reactant  energy. 

When  zero-point  corrected,  this  complex  is  stable  by  only  855  cm*1 .  The  complex  has 
a  triangular  C2v  equilibrium  structure  in  which  the  H-H  internuclear  distance  is  only 
slightly  distorted  from  that  of  H2;  the  B-H  distance  (ca.  2.6  A)  is  much  longer  than  the 
covalent  bond  length  in  BH+  (1.2  A).  Further  along  the  reaction  coordinate,  one  finds  a 
'straight  and  narrow’  reaction  path  characterized  by  stronger  and  stronger  B+-to-H2 
inter-reactant  forces  but  very  small  distortion  of  the  H-H  bond  length  or  the  H-H 
vibrational  frequency.  The  very  restricted  range  of  geometries  (i.e.,  the  narrowness)  of 
this  entrance  channel  and  its  lack  of  ’curvature'  coupling  B+-to-H2  translational 
energy  to  H-H  vibrational  energy  would  be  expected  to  produce  clear  signatures  in  the 
B+  +  H2  ==>  BH+  +  H  experimental  data. 

b.  As  the  B+  ion  approaches  the  H2  molecule  from  very  long  range  (e.g.,  R=  25  A  or 
further),  the  charge-quadupole  interaction  (which  varies  as  R*3)  favors  the  ’T-shaped’ 
C2v  approach  which  eventually  produces  the  above  straight  and  narrow  reaction  path. 
The  permanent  quadupole  moment  of  H2  attracts  positive  ions  to  the  internuclear  regions 
and  repels  positive  ions  from  positions  along  the  H-H  axis  at  long  range.  As  one  moves  to 
smaller  R  values,  the  charge-induced  dipole  interaction  (which  varies  as  R*4)  comes 
into  play.  This  factor  favors  approach  of  a  positive  (or  negative)  ion  along  the  H-H  axis 
(because  a|j  =  0.934  A3  >  a±  =  0.718  A3).  Although  these  electrostatic  and  induced 

interactions  are  dominant  at  very  long  range,  we  find  that  by  the  time  the  B+  is  close 
enough  to  the  H2  moleucle  to  experience  interaction  energies  of  the  order  of  a  few 
kcal/mol,  the  energetically  favored  approach  corresponds  to  a  T-shaped  C2v  structure. 
Although  a  colinear  approach  path  may  have  been  expected  to  be  more  favorable  for 
larger  R  values,  we  find  that  the  linear  structure  becomes  unstable  with  respect  to 
bending  distortion  even  at  rather  large  R  values.  As  a  result,  the  reaction  path  ’bends’ 
toward  the  C2v  path  which  it  then  follows  throughout  the  remainder  of  the  B+  to  H2 
approach. 


c.  The  B+  +  H2  ==>  BH+  +  H  reaction  is  predicted  to  be  endoergic  by  21,318cm'1 
(20,497  cm'1  when  zero-point  corrected).  In  the  guided  ion  beam  experiments33,  no 
flux  of  BH+  product  ions  are  detected  when  B+  ions  collide  with  H2  with  kinetic  energies 
at  or  slightly  above  this  threshold;  this  has  been  used  to  infer  that  an  additional 
activation  barrier  is  present. 

d.  As  one  proceeds  along  the  straight  and  narrow  entrance  channel  starting  from  the  B+ 

+  H2  species  (along  which  C2v  symmetry  obtains  even  though  not  enforced),  a  geometry 
is  reached  at  which  the  asymmetric  distortion  mode  (of  b2  symmetry)  develops  negative 
curvature.  This  occurs  at  an  energy  of 

-25.341096  Hartree,  23,100  cm'1  above  the  B+  +  H2  reactants;  the  geometry  where 
this  occurs  has  an  HH  distance  of  1 .305  A  and  a  BH  distance  of  1 .481  A.  From  here, 
distortion  along  the  asymmetric  mode  leads  to  a  transition  state  that  lies  below  the  C2v  - 
constrained  transition  state.  The  resulting  Cs  transition  state  structure  posesses  one 
'long'  B-H  distance  (rBH=  1 .765  A),  one  'short'  B-H  distance  (rBH'  =  1  -251  A)  and  an 
essentially  'broken'  H-H  bond  (rHH  =  1516  A).  Its  energy  is  23,518  cm*1  above  the 
B+  +  H2  asymptote;  when  zero-point  corrected,  this  point  lies  22,528  cm'1  above  B+ 
+  H2.  This  critical  point  lies  2031  cm'1  or  ca.  0.25  eV  above  the  thermodynamic 
reaction  threshold.  This  energy  gap  corresponds  approximately  to  where  the  guided  ion 
beam  experiments  first  detect  production  of  BH+  product  ions. 

e.  A  geometrically  stable  linear  HBH+^Z)  species  is  found  to  lie  20,892  cm'1 
(19,259  cm'1  when  zero-point  corrected)  below  B+  +  H2.  This  species  lies  on  the 
ground  state  ^1  potential  energy  surface  and  correlates  with  the  B+^S)  +  H2 
reactants  when  C2v  symmetry  is  enforced.  Its  BH  bond  length  (rBH=1-1T3  A)  is  only 
slightly  longer  than  that  in  BH+  (1 .2  A) . 

f.  In  the  absence  of  enforced  C2v  symmetry,  the  1  Ai  and  1 B2  surfaces  of  Fig.  2  are  of 
the  same  (1  A')  symmetry,  and  can  therefore  couple  to  produce  'avoided  crossings'.  As  a 
result  of  these  interactions,  the  ground-state  B+(1S)  +  H2  reactants  can  now  correlate 
directly  to  the  BH+(2I)  +  H  products  as  summarized  in  the  Cs  -  symmety  correlation 
diagram  of  Fig.  4. 


IV.  Discussion  of  Results  and  Summary 

Many  of  the  features  observed  for  the  B+  +  H2  surface(s)  are  remarkably 
similar  to  those  found  earlier  for  Be  +  H2.  In  Fig.  5  we  summarize  the  relative  energies 
of  the  reactants,  transition  states,  and  products  that  arise  in  these  two  systems;  in  both 
cases,  all  energies  are  defined  relative  to  the  ground  state  of  the  M  +  H2  reactants. 

The  energy  of  the  BeH2#  transition  state  corresponds  to  C2v  symmetry  although 
Cs  symmetry  was  explored  in  these  earlier  calculations.  As  Be  approaches  H2,  the 
reaction  path  preserves  C2v  symmetry  and  leads  to  the  transition  state  shown  in  Fig.  5, 
at  which  the  curvature  along  the  asymmetric  b2  distortion  mode  is  positive.  Past  this 
transition  state,  along  the  way  to  the  linear  HBeH  geometry,  a  state  of 1 B2  symmetry 
crosses  the  surface.  At  this  point,  the  curvature  along  the  b2  mode  is  negative.  As  a 


result,  the  reaction  path  moves  away  from  C2v  symmetry  and  produces  BeH  +  H 
products. 

It  may  be  somewhat  surprising  that,  although  the  3P  and  1 P  excited  states  of  B+ 
lie  considerably  higher  than  the  corresponding  states  of  Be  (because  of  the  higher 
nuclear  charge  of  B+),  the  relative  energies  of  the  transition  states,  MH  +  H,  and  HMH 
species  are  rather  similar.  It  should  be  noted  that,  because  the  1  -3P  states  of  B+  lie 
considerably  higher  than  those  of  Be,  the  intersections  of  the  resultant  1  >3B2  states  with 
the  1  A-|  ground  state  are  somewhat  different  for  B+  and  Be.  In  particular,  both  of  the 
1>3B2  states  of  BeH2  intersect  the  1  Ai  state  near  the  C2v  transition  state.  However,  it 
appears  that  the  1B2  state  of  BH2+  does  not  intersect  the  1  Ai  state  as  one  proceeds  along 
the  reaction  path  toward  the  transition  state  (although  it  certainly  couples  strongly  to  it 
to  produce  the  negative  curvature  along  the  b2  mode). 

To  explore  such  differences  and  similarities  further,  it  is  interesting  to 
speculate  about  another  isoelectronic  system:  Li'  +  H2  ==>  LiH~  +  H,  HUH'.  Using  the 
known9  electron  affinities  of  Li  and  of  LiH,  as  well  as  the  LiH  and  H2  bond  energies,  one 
finds  that  LiH'  +  H  should  lie  19,100  cm'1  above  Li'  +  H2.  This  energy  difference  is 
remarkably  close  to  those  for  B+  +  H2  ==>  BH+  +  H  and  Be  +  H2  ==>  BeH  +  H.  The 
energies  of  the  HLiH'  species  and  of  the  C2v  or  Cs  transition  states  for  the  Li'  +  H2 
reaction  are  not  yet  known,  so  further  comparisons  can  not  be  made. 

Major  differences  between  the  Li‘  case  and  those  for  B+  and  Be  involve  the 
location  of  th^  3*1  P  excited  states.  For  B+  and  Be,  the  lowest  of  these  states,  the  3P 
state,  lies  37,300  cm'1  and  22,000  cm1,  respectively,  above  the  ^  ground  state. 

For  Li',  the  3P  state  lies  higher  in  energy  than  Li  (plus  a  free  electron)  and  hence  is 
metastable  with  respect  to  autodetachment.  In  particular,  the  3P  and  1 P  states  are 
believed  to  lie  between  5000  and  10,000 

cm*1  and  between  12,000  and  13,000  cm'1  above  the  state,  while  the  Li'  electron 
detachment  energy  is  only  5000  cm'1.  For  these  reasons,  the  intersection  of  the 
ground  state  C2v  surface  for  Li"  +  H2  by  the  excited  1 B2  surface  is  expected  to  occur  at 
much  lower  energy  than  in  the  B+  and  Be  cases.  Moreover,  the  2Ai  surface 
corresponding  to  the  autodetached  Li  +  H2  species  is  also  expected  to  come  into  play  at 
low  collision  energies.  Therefore,  qualitatively  different  behavior  is  expected  both  in  the 
low-energy  potential  surfaces  of  Li'  +  H2  and  in  the  guided  ion  beam  experiments 
studying  Li'  collisions  with  H2  .  It  is  our  intention  to  explore  this  interesting  case  in 
the  near  future. 

Before  closing,  it  should  be  mentioned  that  there  are  interesting  aspects  of  the 
B+  +  H2  potential  energy  surfaces  that  were  not  addressed  here  because  of  our  emphasis 
on  the  lowest  singlet-state  potential.  In  particular,  the  location  of  the  3B2  surface  in  the 
neighborhood  of  the  singlet  state's  Cs  transition  state  is  of  substantial  importance  to  a 
full  interpretation  of  the  B+  +  H2  reactivity  data.  If  transitions  to  the  1 B2  surface  play 
an  important  role,  transitions  to  the  corresponding  triplet  surface  will  also  be 
operative.  The  latter  events  may  occur  with  reduced  probability  because  of  the  need  for 
singlet-triplet  coupling,  but  they  will  occur  at  lower  energy  because  the  3B2  state  lies 
below  the  1 B2  state. 

Finally,  in  this  work  emphasis  has  been  placed  on  the  reaction  path  as  it  enters 
the  'straight  and  narrow’  channel  of  C2v  symmetry;  little  has  been  said  of  colinear  M- 
HH  geometries.  Hurst  has  shown  that  the  linear  approach  does  not  lead  to  a  lower-energy 


path  to  the  formation  of  BH+(2i)  +  H,  and  our  data  confirms  that  the  energetically 
favored  approach  evolves  into  the  C2v  geometry  once  the  B+  and  H2  begin  to  interact  via 
chemical  valence  forces. 
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Table  I.  Total  and  Relative  Energies,  Geometries,  and  Vibrational  Frequencies  for  Species 
Along  the  B+  +  H2  ==>  BH+  +  H  Reaction  Path3 


Species 

Method 

Used 

Electronic 

Energy 

(Hartrees) 

Optimized 
Geometry  (A) 

MCSCF 

Frequencies0/ 
Zero  Point 
Energies 
(cm*1 ) 

Relative 

Energies 

(cnr1)b 

B+^S)  +  H2 

MCSCF 

CCSD(T) 

-25.446357 

-25.468865 

rHH  =  0.755 

OHH  =  0.734 

4224/2112 

0 

bh2+  C2v 

complex 

MCSCF 

-25.468239 

rHH  =  0.762; 
reH=  2.605 

230  (ai), 

437  (b2), 

CCSD(T) 

-25.474069 

rHH  =  0.743; 
rBH=  2.605 

4133  (ai )/ 
2400 

-1143 

bh2+  C2v 
transition 

mcs  ;f 

-25.322627 

rHH  =  1-396; 
rBH=  1.411 

451 2i  (aO, 
1279  (ai ), 

state 

CCSD(T) 

-25.352230 

3424i  (b2)d 

25,595 

HBH+(1I) 

MCSCF 

CCSD(T) 

-25.524859 

-25.563931 

rBH=  1.183 
rBH‘  =  1-172 

2911  (ai), 
2632  (b2), 

973 

(bending)/ 

CCSD(T) 

-25.564067 

rBH=1  -173 

3745 

-20,892 

bh2+cs 

transition 

state 

-25.332540 

-25.361695 

rHH  =  1.516 
rsH=  1-765 
tbh’  =  1-251 

2840i, small®, 
2087/  ca. 

1050 

23.518 

BH+(2x)  + 

MCSCF 

-25.351373 

rBH  =  1.199 

2582/ 

H 

CCSD(T) 

-25.371722 

rBH  =  1-198 

1291 

21,318 

BHH+  linear 
complex 

MCSCF 

CCSD(T) 

-25.446663 

-25.469942 

rHH  =  0.756 
rsH=  2.729 

112,  4202, 
303i 

(bending)d 

CCSD(T) 

-25.488630 

rHH  =  0.702 
rBH  =  2.865 

/2 1  57 

-4337 

a.  Where  separate  geometry  optimizations  were  carried  out  at  the  MCSCF  and  CCSD(T) 
levels,  two  sets  of  geometries  are  reported.  In  all  cases,  the  MCSCF  geometry  is  listed 
first,  and  the  CCSD(T)  geometry  appears  second.  Where  geometry  optimization  could  not 
be  carried  out  at  the  CCSD(T)  level,  the  CCSD{T)  energies  were  computed  at  the  MCSCF 
geometries. 

b.  In  all  cases,  only  the  CCSD(T)  energies  are  used  because  they  represent  our  best 
values.  They  are  given  relative  to  the  B+  +  H2  reactants.  These  are  electronic  energies, 
and  thus  do  not  include  zero-point  corrections. 

c.  These  local  harmonic  frequencies  were  obtained  from  the  analytical  second  derivatives 
of  the  MCSCF  energy  at  the  MCSCF  geometries. 

d.  It  is  not  appropriate  to  compute  zero  point  energies  here  because  this  is  neither  a  true 
minimum  nor  transition  state  (i.e.,  more  than  one  imaginary  frequency  appears). 

e.  The  precise  location  of  the  Cs  transition  state  was  very  difficult  to  determine.  See  text 
for  further  discussion. 


Figure  Captions 


Figure  1.  Two-dimensional  contour  plot  of  the  ground-state  1  Ai  potential  energy 
surface  for  the  Czv  insertion  of  Be{1S)  into  H2  to  produce  linear  HBeH(1Z).  The 

contour  spacings  represent  approximately  5,000  cm'1  in  energy.  The  labels  a-d,  f 
refer  to  the  geometrical  points  discussed  in  Sec.  I.C.  Along  the  vertical  axis  is  the 
distance  from  the  B+  to  the  center  of  the  H-H  bond;  the  horizontal  axis  labels  the  H-H 
distance. 

Figure  2.  Configuration  correlation  diagram  for  C2v  insertion  of  B+(1S)  into  H2  to 
produce  linear  HBH+.  The  energies  are  in  untis  of  1000  cm'1.  The  lai,  2ai,  and  1b2 
orbitals  correspond  to  the  H2  og,  B+  2s,  and  B+  2p  (in  plane)  orbitals,  respectively  for 
the  B+  +  H2  reactants.  For  the  HBH+  product,  the  lai  and  fb2  orbitals  are  the  two 
(symmetric  and  assymmetric)  c  B-H  bonding  orbitals. 

Figure  3.  Constructive  interaction  between  the  in-plane  M  np  orbital  of  b2  symmetry 
and  the  antibonding  H2  ou  orbital  also  of  b2  symmetry. 

Figure  4.  Configuration  correlation  diagram  for  Cs  insertion  of  B+(1S)  into  H2  to 
produce  either  linear  HBH+  or  BH+(2S)  +  H.  The  energies  are  in  units  of 
1000  cm'1.  The  la-,  2a*.  and  3a'  orbitals  correspond  to  the  H2  og,  B+  2s,  and  B+  2p 
(in  plane)  orbitals,  respectively  for  the  B+  +  H2  reactants.  For  the  HBH+  product,  the 
1  a*  and  3a’  orbitals  are  the  two  (symmetric  and  assymmetric)  o  B-H  bonding  orbitals. 
For  BH+  +  H,  la'  is  the  BH+  a  bonding  orbital,  3a’  is  the  BH+  non-bonding  c  orbital, 
and  2a'  is  the  H  atom  Is  orbital. 

Figure  5.  Energy  diagram  showing  relative  energies  of  M  +  H2,  HMH,  MH  +  H,  and 
transition  states  for  the  B+  +  H2  and  Be  +  H2  cases.  Also  shown  are  the  locations  of  the 
1'2P  excited  states  of  the  B+  and  Be  species.  All  energies  are  in  untis  of  1000  cm'1. 
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